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Abstract 



The derivative of the associated Legendre function of the first kind of integer degree with 
respect to its order, dP£(z) /dfi, is studied. After deriving and investigating general formulas 
for /i arbitrary complex, a detailed discussion of [dP% (z)/dfj]fj,=± m , where m is a non-negative 
integer, is carried out. The results are applied to obtain several explicit expressions for the 
associated Legendre function of the second kind of integer degree and order, Qn m {z). In 
■ particular, we arrive at formulas which generalize to the case of Qt m (z) (0 ^ ra n) the well- 

' known Christoffel's representation of the Legendre function of the second kind, Q n (z). The 

Q\ , derivatives [d 2 P£(z)/dfi 2 ]^ =m , [dQk(z) / '9/x] M=m and [dQt n _ 1 (z) /d/i] M=m , all with m > n, 

&\ ' are also evaluated. 



o 
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1 Introduction 



KEY WORDS: Legendre functions; parameter derivative; special functions 
AMS subject classification: 33C45, 33C05 



It is the purpose of the present paper to contribute to the theory of special functions of mathe- 
matical physics and chemistry. Specifically, we shall add to the knowledge about the associated 
Legendre functions of the first, P^{z) 1 and second, Q£(z), kinds (cf, e.g., [1-29]). We shall touch 
two particular problems. First, we shall investigate the derivative of the associated Legendre func- 
tion of the first kind of integer degree with respect to its order. Second, it will be shown that the 
results of that investigation may be used to construct, in a straightforward and unified manner, 
several known representations of the associated Legendre functions of the second kind of integer 
degree and order, Q™(z); in the past those representations were obtained by other authors with 
the use of a variety of, usually more complicated, techniques. In addition, some possibly new 
expressions for Q™{z) will be also presented. 

The literature concerning the derivative dP^(z)/dfi is very limited. Surprisingly, no relevant 
expressions have been given in [29] , which otherwise contains a large collection of parameter deriva- 
tives of various special functions. Several variants of the formula 

dPftz) 



dji 

where 



= iP(v + l)P„(z) + Q v (z), (1.1) 



1 



is the digamma function [14,20,23,25,30], may be found in [23, page 178]. Robin [18, Eq. (333) 
on page 175] gave the following representation 1 of dP^{z)jd\x: 



d/j, 2 v w z - 1 

k\T(v - k + l)T(k - n + 1) 



fc=0 



(1.3) 



Brown [31], being apparently unaware of the Robin's finding, rederived the expression in Eq. (1.3) 
in the particular case of v = n G N and used it to find a representation of Q™(z) suitable for 
numerical purposes. An application of the derivative dPj?(z)/d(i, to the evaluation of the integral 

J dz [P^ (z)] 2 1 (z 2 — 1), where C is the contour starting from z = and returning to it after a 

positive circuit round the point z — 1, was presented by Watson [32]. 

The detailed plan of the paper is as follows. In section 2, we shall summarize these facts 
about the associated Legendre functions of the first and second kinds of integer degrees which 
will find applications in later parts of the work. In section 3, the derivative of the associated 
Legendre function of the first kind, dPfi(z)/dfj, (with z G C \ [—1,1]), is investigated. After 
general considerations forming section 3.1, where no restrictions are imposed on the order \x of 
the Legendre function, in sections 3.2 and 3.3 we focus on the derivatives [dP^(z)/dfj]^ = ± m with 
m G N. The cases of ^ to ^ n and m > n arc treated separately A brief discussion of 
the derivative dP£(x) / 'd\x (with x G [—1, 1]) is contained in section 3.4. Section 4 presents some 
applications of the results of section 3. In section 4.1, we obtain several explicit representations 
of the associated Legendre function of the second kind of integer degree and order, Q^ m (z) and 
Qn" l ( x )- Some of these representations seem to be new, particularly these ones which generalize 
the well-known Christoffel's formulas for Q n (z) and Q n (x). Finally, in sections 4.2 and 4.3, the 
derivatives [d 2 P^(z)/dn 2 } ll=mi [dQ^(z)/d/j]^ =m , and [dQt n _ 1 (z)/dfi] ti=m , all with m > n, are 
expressed in terms of the derivatives [<9P^(±z)/ctyi] M= _ m . 

Throughout the paper, we shall be adopting the notational convention according to which 
v G C, \i G C, n G N, m G N, x G [-1, 1], and z G C \ [-1, 1], the latter with 

- 7r < arg(z) < 7r, — 7r < arg(z ± 1) < 7r, (1.4) 

hence, 

-z = e* iw z, -z + l = e T ' m (z-l), -z - 1 = e T ' m {z + 1) (arg(z) ^ 0). (1.5) 

Furthermore, it will be understood that if the upper limit of a sum is less by unity than the lower 
one, then the sum vanishes identically. 

The definitions of the associated Legendre functions of the first and second kinds used in the 
paper are those of Hobson [12]. 



2 The associated Legendre functions of integer degrees 

The material presented in this section, to be referred to in later parts of the paper, is based 
primarily on [12, 17, 18]. 



1 The notation used in the present work differs in some respects from that adopted by Robin in [17-19]. In 
particular, the digamma function used by Robin was defined as 

^(0= 1 dr(c + 1) 

r(c + i) dc 

rather than as in our Eq. (1.2). Also, one should be warned that in Eq. (333) on page 175 in [18] the factor p\ in 
the denominator of a summand is missing. 
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2.1 The associated Legendre function of the first kind of integer degree 

The associated Legendre function of the first kind of a non-negative integer degree n and an 
arbitrary complex order /i may defined as 



or equivalently as 

D „, . 1 T{nTV+l) (z + iy /2 (z±lY' ( z T l 



r(iT M )r(n- M +i) v 2 / V '^'«±i 



(2.2) 



with either upper or lower signs chosen. Both hypergeometric functions in Eq. (2.1) are proportional 
to the Jacobi polynomial [23,25] P^T >J '' IJ '\z); in terms of the latter, one has 

Combining this with the Rodrigues-type formula for the Jacobi polynomials, which is 

1 H n 

P ^ )(Z) = 2^ {Z iyU{z + ^dz^ [(Z ~ l)n+a{z + ir+/3] ' (2 ' 4) 

leads to the following Rodrigues-type representation of P£{z): 

A*/2 



*<*> - 2T ( „', + i) (St) £ K' - »t + "-I 



(2.5) 



The function P%{z) is single- valued in the complex z-plane with a cut along the real axis from 
— oo to 1. For a negative integer degree, the associated Legendre function of the first kind is defined 
through the relationship 

Ptt n -M = P£{*)- (2-6) 
For the sake of later applications, it is convenient to rewrite Eq. (2.1) as 



and Eq. (2.2) as 

P£(z) = n!r(n + /i+l) 



z + lY /2 (z±\ 



z-l 



^ k\(n - k)\T(k T M + l)r(n - k ± fi + I) \z±l' ' ^ 



fe=0 



From either of Eqs. (2.1), (2.2), (2.7) or (2.8), the function P£(z) is seen to possess the reflection 
property 

= ^Y^TTf^- < 2 - 9 > 

If fi = ±m, then it holds that 

P™(z)=0 (m>n), (2.10) 
P r f m (-z) = (-TP± m (z) (0 < m < n), (2.11) 
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p - m(z) = {ih^f Fr(z) (° ( 2 - 12 ) 

If < m n, from Eqs. (2.7) and (2.8) wc have 

^ 4 J fc!(fc + m )!(„- m -jfc)! ^ 2 ) {U ^ m nh 

(2.13) 

pm(A , , n (n + m)\ fz + l\ m/2 " fc (fc + n)! + ^ ; 

" 1 ; " k J (n-m)!V«-iy £^ ; k\(k + m)\{n-k)\ \ 2 

(0<m<n), (2.14) 



P™(z) = n!(n + m)! 



z=Fl\ 7 /z ± 1 



z± 1 / V 2 



n — m ^ ✓ ^ ^ s fc 

*E fe!(fc + m)!(n-fe)!(n-m-ifc)! (Slj <° < m < »)' ( 2 ' 15 ) 

and also 

/ x | / 2 -i \ m / 2 

P-™(Z) = (-)n+mt n - m )- ( Z 



(n + m)\ \ 4 

*E (-)* ,,,, (fc + " + m)! M , (0<m<n), (2.17) 

^ v ; fc!(fc + m)!(ra-m-fc)! V 2 / 



fc=0 



,-<„ . „ !( „_ m)! (-i)"' 2 (i|i)" 

1 

fc!(fc + m)!(n — fc)!(n — m — fc)! ± 1 



n — m _^ / 1 \ ^ 

X y 7777 777 T777 777 \ ~~7~T ) (° < TO < «)■ ( 2 - 18 ) 

fe=0 



For m > n, the same equations yield 



™ - (^l)^t M(t ^ff_ t) , (i^)' <„>„>, (, 19 , 



P " ^ (n + m)!(m-n-l)! U + £ fc!(n-fc)! V 2 

(to > n), (2.20) 

n! /z-l\ m/2 /z+l\ n A fc (fc + TO-n-1)! /z-l\ fe 

" lJ " (m-n-1)! \z + l) V 2 / ^ + m)!(n - fc)! \z + 1) 



(m > n), (2.21) 
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and 



Pn" 1 ^) = (-)" 



z - 1 



^/2 



(m — n — 1)! \0 + 1, 
fc (m-fe-1)! 



z - 1 



EH 

fe=0 



fc!(n — fc)!(ra + m — fc)! \z— 1 



z + 1 



(m > n). (2.22) 



In deriving Eqs. (2.20) to (2.22), the formula 

lim + =( _)nW (m >n,n') (2.23) 

n—>-m r(n + /x + 1) (m-n-1)! 

appears to be useful. 

Following Hobson [12], on that part of the cut for which z = x(— 1 ^ a; ^ 1) we define the 
associated Legendre function of the first kind of non-negative degree as 

P£{x) = i [e i7r ^ 2 P£(x + 10) + e-' ln ^ 2 P£(x - i0)] = c ±i7r ^ 2 P£(x ± iO). (2.24) 

2.2 The associated Legendre function of the second kind of integer de- 
gree 

The associated Legendre function of the second kind of a non-negative integer degree n and an 
arbitrary complex order fi may be defined by either of the two expressions: 



it e' 



2 sin(7r/x) 



_ i> + m + i) , 

" l j r(n- M +l) " 1 ' 



(2.25) 



(2.26) 



2 sin(7r^) 

equivalence of which follows immediately from Eq. (2.9). An extension of the definition to negative 
integer degrees is made via the formula 



Q- n -i(z) = -Q*{z) + n-rj-rPXiz). 



One obtains 



and 



it e 



2 sin(7r/i) 



sin(7r^) " 
T(n - jji + 1) 



Q*„_i(*) = 



7r e 



2 sin(7r^) 



KW + (-) n ^M]. 



From Eq. (2.25) it is seen that 



Q n W e r(n + ((+1) W 



(2.27) 

(2.28) 
(2.29) 

(2.30) 



It should be observed that the function Qt n _ 1 (z) does not exist if yU — n is a negative integer or 
zero. The functions Q£(z) and Q^„_ 1 (z) (provided the latter exists, see above) are single-valued 
in the complex plane with a cut along the real axis from — oo to 1. 
For ^ m ^ n, it is known 2 (see [18, pages 81, 82, and 85]) that 



Qn( z ) = \Pn{z) In i±l - W^-iC*) (0 < m < n), 



(2.31) 



2 In [18, pages 81, 82, and 85] expressions for Q™(z), equivalent to these following from our Eqs. (2.31) to (2.33), 
have been given in somewhat different forms (in this connection, cf footnote 1 on page 2). We have modified Robin's 
formulas to make them concurrent with the notation used in the present paper. 
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with 



C-iW - ^(n + m+l)+V>(n-m + l)]fTO 



(-) m fz + l\ m/2 "^^ Nfc (fc + n)!(m-fc-l)! 



z-17 ^ v 7 fc!(ra-fc)! 

fc=0 



l(n + m)! ^z- l^ m/2 ^ (fc + n)h/>(fc + m + l) 



2(n-m)! V-2 + 1/ ^ fe!(fc + m)!(n - fe)! \ 2 

™/2"-m (fc + n + m )!^(fc + l) /z-l\ fe 

g o fc!(fc + m)!(n-m-fe)! (—J (° < ™ < ») 



(2.32) 



or 



ICiW = ~[V(n + m+l) + V(n-m+l)]^) 



(-)"+" /z-ir^y, 1 (fc + n)!(m-fc-l)! fz + l x fe 



^ fe!(n-fc)! 

fc=0 



(-)" (n + m)! (z + l\ m/2 fc (fc + n)h/>(fc + m + l) /z + lV 

2 (n-m)!V«-l/ ^ M fc!(jfe + m)!(n - fc)! V 2 / 

R"+ m ^z 2 -lV" /2? ^, ^ fc (fc + n + m)h/-(fc + l) ^z + l^ fe 



\ ffl/ £ ft — ffl 

■) Eh 

/ z — n 



fc!(fc + m)!(n — m — k)\ \ 2 
(O^m^ n). (2.33) 

Another expression for W^ l _ 1 (z) may be deduced from the findings of Brown [31]; it is 
{-) m (z + l\ m/2r ^ k (k + n)\{m-k-l)\ fz-1 



w »-m - ~—\—) £ H — w=w — {- 

(-) n+m /z-ir' 2 ^ 1 (k + n)\(m-k-l)\ (z+\\ k 

+ 9 I v 4- 1 i 



2 V^ + iy ^ fc!(n-fe)! V 2 

fc=0 



fc!(fc + m)!(ra-m-fc)! V 2 



R" +m / z 2 -l V" /2T y^ fc (fc + n + m)!^(fc + l) f z + 1 
+ 2 V 4 / ^ fc!(fc + m)!(n-m-fc)! V 2 
(0<m^ n). (2.34) 

Furthermore, Snow [13, pages 55 and 56] gave a representation of W^ l _ 1 (z) which may be easily 
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shown to be equivalent to 

1 ( z =p l\ m/2 ( z ± i\""~ m 

W ^ {Z) = ± 2 nl{n + m)l {—l) {—) g fcl(fc + m)l(n-fc)l(n-m-fc)l 



1 



x [V>(n -m-k + l) + ip(n - k + 1) - ip(h + m + 1) - + 1)] 



±-n!(n + m)! , 
2 V^T 1 



1 1 \ m/2 / 1 

2 ± 1 \ / z =F 1 



fc=i 



(fe-l)! 



£Tl 

(fc + n)\(k + n- m)\(m - k)\ \z±l 



=F ' n!(n + m)! 



2 

m— 1 



l/2 



Z± 1 



\- / )k (m-k-l)\ ( zjfl 

' k\(n-k)\(n + m-k)\\z±l 

k—0 x 



(0 < m < n). 



(2.35) 



Exploiting Eq. (2.31) and either of Eqs. (2.32) to (2.34), Q n m (z) may be evaluated using 

(»-"0»™^ ,n^_^ (2 . 3(;i 



(n + m)! 



Q™(z) (0 < m < n), 



which is the direct consequence of Eq. (2.30). The functions Q z t.™_ 1 (z) with ^ m ^ n do not 
exist. 

For m > n, it holds that [18, Eq. (63) on page 35] 



Quiz) 



(-)' 



-(n + m)!(m - n - 1)! [P„- m (-z) - (-)"p-" l (z)] (m > n). 



(2.37) 



The function Q~ m (z) with m > n does not exist. 

On the cut x G [—1,1], again following Hobson [12], we define the associated Legendre function 
of the second kind of non-negative and negative integer degrees as 



and 



respectively. 



e -i*M/2Q£(a; + io) + e^/ 2 Q^{x - iO) 
'e-^Qtn-dx + iO) + e i7r ^ 2 Q^„_ 1 (a; - iO) 



(2.38) 
(2.39) 



3 The derivatives dP£(z)/dfi and dP%(x)/dii 
3.1 General considerations 

We begin with the observation that after differentiating Eqs. (2.6) and (2.9) with respect to fi, one 
arrives at the following two general relations involving the derivative in question: 



dPH n _M _ dP%{z) 



(3-1) 



-hHn + M + 1) + V(n — M + + (-) 



n r(n- M +l)^(-z) 
r(n + /x + 1) 5/i 



(3.2) 



7 



Because of Eq. (3.1), in the remainder of the paper we shall focus on evaluation of dP£(z)/d/j,. 
We proceed to the construction of several explicit representations of dP^(z)/dfi. 
From Eq. (2.5), we have the Rodrigues-type formula 



= -^WlnJij+V(n-M+l)^W 



+ - 



1 (z - 1 

2 n T(n-n + l) \z + l 

This may be rewritten as 



M/2 



d n 

dz™ 



(z- l) n -»{z + 1)"+" In 



z + 1 
z-1 



(3.3) 



(3.4) 



with 



-P;(z)ln^+^^ + l)TO 



+ 



2 n T(n- fj,+ l) V-z + 1 



A</2 



dz n 



{z-l) n -»{z + 1)"+" In 



+ 1 



z-1 



(3.5) 



Choosing in Eq. (2.7) the upper set of signs and differentiating with respect to yields 
dP%(z)/dn in the form (3.4) with 



+ 1 y/2 n ( fc + n )|^( fc _ M+1 ) / z _l 



(z+iy^ (fc+7 



k)T{k-n + l) V 2 



(3.6) 



An alternative representation of U%(z) is obtained if in Eq. (2.7) one chooses the lower set of signs 
and then differentiates with respect to ^; this results in 



UH(z) = [^(„ + ^ + i)+V>(n-/i+l)]P£(*)-(-y 



,r(n + /i + l) /z + 1 



r(n-/x+l) V-2- 1 



p/2 



Eh 

fe=0 



fc (fc + n)h/>(fc + pt+ 1) ^z + 1 
fc!(n - k)\T(k + fi+l) 



(3.7) 



Playing with Eq. (3.7) with the aid of Eq. (2.7) (with the lower signs chosen) and of the following 
easy to prove identities: 



£(n + M+l) = +n , r(-n'- M ) 



T(n' + /i+l) v y r(-n-/i) 
V>(n + //. + 1) - -0(n' + /i + 1) = ^(-n - m) _ V"(- n ' - M) 



and 
leads to 

Ltf(z) = [^(n- M +l)+V<-n-M)]TO- 



(3.8) 
(3.9) 



z + 1 

r(n - /i + i)r(-n - /i) V-z - 1 



H/2 



E 

fc=0 



(k + n)!r(-fc - - AO A + 1 



fc!(n-fc) 



(3.10) 



Furthermore, differentiating Eq. (2.8) with respect to /i gives 



- 1 p(n + l x+l)PH(z)±nT(n + f x+l) 



z + 1 
z-1 



M/2 



Z± 1 



ip(k T V + 1) ~ ip(n - k ± fi + 1) (z^l 
^ k\(n - k)W(k + m + l)r(n - fc ± /i + 1) \z± 1 



E 



(3.11) 
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Manipulations with Eq. (3.11), similar to those which have led us from Eq. (3.7) to Eq. (3.10), give 



r(— n — (j,) \z — 1 



z+1 



z + 1 



EH* 



T(k — n — fx) 



k=0 



fc !(n-W-M + l) ^ (fc " M+1) "^" n " M)] l^TT 



z - 1 



(3.12) 



or 



n 

><E(-) fe 



r(-n - m) 

r(-fc- M ) 



M/2 



Z - 1 



fc=0 



ft!(ra-fc)!r(n-fc-/z + 



-[^(n-fc-M+^-V'C-fc-M)] 



z - 1 



(3.13) 



according to whether the upper or the lower signs are chosen in Eq. (3.11). 

Next, let us look at Eqs. (2.7) and (3.6) which give the functions P£(z) and U%(z), respectively. 
It is seen that both these functions are of the form [(z + l)/(z — l)]^/ 2 times a polynomial in z — 1 
of degree n [except for the case of \x — m > n, when P™(z) vanishes identically, see Eq. (2.10)]. 
Consequently, if /i ^ m, it must be possible to represent U%(z) in the form of a linear combination 
of the functions Pj^(z) with degrees not exceeding n: 



(3.14) 



fc=0 



[in the case of /i = m one may seek U£{z) in the form of an expansion analogous to (3.14) but with 
Pg(z) replaced by P^(-z); cf alsoEq. (3.34)]. Comparing coefficients at (z-l)"[(z+l)/(z-l)]^/ 2 
in Eqs. (2.7) and (3.6), we see that the coefficient c^„ in the above expansion is 



ij){n — [i + 1). 



(3.15) 



The reasoning leading to determination of the coefficients c^ k with 0^/c^n— lisas follows. 
If we differentiate the well-known Legendre identity [23, 25] 



_(l-, 2 )_ +n (n + l) 



1-z 2 



PH(z) = 



with respect to /u, this yields 



d /, 9n d / , n 
d- Z {1 - z) d- Z +n{n+1) 1-z 2 



dP£{z) 
dfi 



2fi 



On exploiting the identity (3.16), we also obtain 



d ,„ 9s d . „. u 2 



(3.16) 



(3.17) 



(3.18) 



By subtracting Eq. (3.18) from Eq. (3.17), we find that the function U£{z) satisfies the inhomoge- 
neous differential equation 



d d , „, u 2 



UH(z) = -2 



dP%(z) , 2/i 



dz 



+ 



1 -z- 



(3.19) 
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For a while, let us focus on the expression on the right-hand side of Eq. (3.19). We shall show 
how it may be developed into a finite sum of the functions P£(z) with ^ k ^ n — 1. To begin 
with, we observe that it holds that 



dz 



1 - z 



(z-1) 



dz 



z + 1 



PX(z)-nP£(z) 



{z + 1} ^£) + Jt-PH{z) nP^(z) 



If use is made of the well-known relation [23, 25] 



{z 2 1)^1 = nzPZ(z) (n + ^P^z), 



(3.20) 



(3.21) 



the expressions in the square brackets appearing on the right-hand side of Eq. (3.20) may be written 



as 



(*Tl) 



dz +lfl^)-^) = =F 
Consider next the recurrence relation [23,25] 

(2k + l)zP£(z) = (k-n+ l)Plt +1 (z) + (k + riP^z). 
Evidently, it may rewritten as 

(2k + i)(z + i)p£(z) = [(k-» + i)p>; +1 (z) + (k + » + i)pj;(z)} 

+ [(k- /J )P£(z) + (k + i i )P£_ 1 (z)]. 

Multiplying both sides of Eq. (3.24) by (— ) fe+ ™, summing over k from fc = 0toA; = n— 1, utilizing 
then the property [cf Eq. (2.6)] 

Pg(z) = P»_ 1 (z) 

and dividing the result by z + 1 , wc find 



(«T/«(s)±(n + 0)f* , _ 1 (*) 



(3.22) 
(3.23) 



£(-) fc +"(2fc + l)P£(*) 



(n- M )P£(z) + (n + M )P„%(z) 



fe=0 



z + 1 



If on both sides of Eq. (3.26) use is made of Eq. (2.9), the former becomes 



V(2fc I n E(fc+ii+i) p- 
^ 2fc + 1 ) r( *- /1 + i) p * 



'(-*) = - 



r(n + M +l)^(-^)-^i(-^) 



r(n - m) 



z + 1 



From this, after replacing /x with — /i and z with -z, we infer 

n-l 



fe=0 

Hence, it follows that 

n-l 



+ ^rf+T+T) fe w " r(n + M) z - ! • 



r(n-/i+l) 



E( 2fc+1 ) 



fc=0 



r(fc-// + 1) 
r(fc + Ai + 1) 



r(n + it) 

(n + M )P^(z)-(n + M )P^ 1 (z) 
z-1 



(3.24) 
tilizing 

(3.25) 
(3.26) 

(3.27) 

(3.28) 

(3.29) 



Combining Eqs. (3.20), (3.22), (3.28), and (3.29), we arrive at the sought expansion 

2/u 



, dP£(z) + 



dz 



1-z 2 



TO 



= ]T(-) fc +"(2fc + l) 



fc=0 



l-H 



fc+n r(rc + M + i)r(fc -/* + !) 

r(n-it+l)r(fc + M+l) 



to, 



(3.30) 
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which, apart from being useful in the present context, seems to be also interesting for its own sake. 

We return to the problem of determination of the coefficients c^ fc . In virtue of the result in Eq. 
(3.30), we rewrite the differential relation (3.19) as 



d ,„ 9x d , 
d^ (1 -^ +n(n+1) 



n-l 



£(-)*+»(2fc + l) 



k=0 



l-(-) 



fc+n i> + ^ + i)r(fc- M + i) 

r(n - n + l)T(k + n + 1) 



P£{z). (3.31) 



On substituting the expansion (3.14) into the left-hand side of Eq. (3.31), after equating coefficients 
at Pfr(z) on both sides of the resulting identity, we obtain 



C nk 



k+n 



(n - k)(k + n + 1) 
(0 sC k ^ n- 1). 



2/,- + i r _ fc+n r(n + M + i)r(fc-/i + i) 

1 ' T(n - n + l)T(k + n + 1) 



(3.32) 



Plugging Eqs. (3.15) and (3.32) into Eq. (3.14), we arrive at the following representation of U£{z): 

71-1 



fe=0 



2k + 1 



i-H 



(n- fe)(fe + n+l) 

fc+» r ( n + ^ + i)r(fc-M + i) " p ^ 



r(n- /U + i)r(fc + M + i)_ 

With the use of the property (2.9), the above result may be transformed into 



n-l 

+E(-) fe+ " — 

k=0 



2k + 1 



(n - k)(k + n+l) 



fc w v ; r(n-/i+l) * 



The last expression is of value when ^ yU = to ^ n. 

3.2 Evaluation of [9P^(^)/9/i] M=m 
3.2.1 The case of < to < n 

For = m, the fundamental equation (3.4) becomes 



5P£(z) 



fi—m 



= Ip^^ln J±l + (0 < m < »). 



It follows from Eq. (3.5) that the Rodrigues-type representation of U™(z) is 



-JTO In i±| + ^(n - m + l)i^(z) 
z — 1 



+ 



1 



z - 1 



2"(n - to)! V- 2 + 1 
(0 ^ to ^ n). 



m/2 



dz r ' 



(z-l)"- m (z + l) 



n+m 



In 



z- 1 



If in Eq. (3.6) use is made of the well-known limiting property 

lim m~"tn = H fe+ro (™ - fc - !) ! > k ) 

H^m I (k — /I + 1) 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



(3.37) 
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and if, whenever necessary, the gamma function is replaced by the factorial, after some straight- 
forward manipulations we obtain the following representation of U™(z): 



Tjmj y \ , xm f z + l \ m/2 ^ u (fc + n)!(m-fc- 1)! (z - ■! 

ni) ( } U-J to } k\(n-k)\ \ 2 



fz 2 -l\ m/2n ^ (fc + n + ro)ty(* + l) ( z - \\ k 

+ {—) ]2 k Kk + m)K n- m -ky. {—) 

(3.38) 

Similarly, from Eq. (3.7) it follows that 

U™(z) = lyj(n + m +l)+yj(n-m+l)}P™(z) 

(n + m)! /z + l\ m/2 " (fc + n)!^(fc + m + l) /z + l\ fc 
1 J (n-m)!V«-l/ ^ J fc!(jfe + m)!(n-fc)! V 2 / 
(0 < to < n). (3.39) 



Furthermore, from Eq. (3.11) one finds 

EC(*) = 4>{n + m+l)P™{z) + {-) m n\(n + m)\ 

E" i \k (m — fc — 1)! /. 
_ M fc!(n-fc)!(n + m-/c)! V 



-i \ m/2 / -. 
Z + 1 \ ' / Z + 1 



Z - 1 



— n!(n + to)! 



-. \ m/2 / -. 

z — 1 \ / z + 1 



z + 1 



^ k\(k + m)\(n - k)\(n - m - k)\ \z+lj v 7 v ; 



and 

t/™(z) = ^(n + TO + l)P„ m (z)+n!(n + TO)! 



, x m/2 / , -, 

z — 1\ I z + 1 



z + iy V 2 

fc 



"t {k + n)\{k + n-m)\(m-k)\\z-lj 
+nl(n + m)\ 



. , \ m/2 , 1 
z + 1 \ ' / z — 1 

z-l r 

E m ip(n - to — k + 1) — ^(fe + to + 1) /z + 1\ fc 
fc!(fc + TO)!(n-fc)!(n-TO-fc)! ^^ij (U ^ m ^ n). 

(3.41) 

Finally, from either of Eqs. (3.33) or (3.34), with the aid of Eqs. (2.13) to (2.15), one deduces 
that 

ITW - «» - m + Wl - M g H > (n _ » + l - t) „-.(-,) 



(3.42) 
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or equivalcntly 



U?(z) = V(n-m+l)P„ m (z)- 



Eh* 



2k + 1 



— m — 1 

+ E (-) 



(n + m)! 

(n — m)!^^ v ; (n — fc)(fc + n+l) 
2fc + 2m + 1 



fc=0 



(n — m — k){k + n + m + 1) 



i-H 



k+n+m 



k\(n + m)\ 



(k + 2m)!(n - m)! 



P fe ™ m (z) (0<m<n). (3.43) 



3.2.2 A closer look at the case of m = 
Setting n = in Eq. (3.17) results in 



_(!-*»)_ +„(„ + !) 



ap^(z) 



= 0, 



(3.44) 



A»=0 



i.e., [dP%(z)/d[j]n=o solves the same differential equation as the Legendre polynomial P n {z) = 
P®{z) does. Consequently, it must be possible to write it in the form of a linear combination of 
P n (z) and the Legendre function of the second kind Q n (z) = Q°(z). 

To find that combination, at first we observe that for m = Eq. (3.35) becomes 



dP£(z) 



fi=0 



\p n {z) In £±I + tf°( z ). 



(3.45) 



From Eq. (3.43) we see that U®(z) is a polynomial in z of degree n, which may be written as 



U°(z)=iP(n + l)P n (z)-W n - 1 (z). 
Here W n -i(z) is a polynomial in z (of degree n — 1) given by 

1 / \fc+n 



1 - C-r 

W n . 1 (z) = 2^ V" 



2k + 1 



or equivalcntly 



fc=0 



int[(n-l)/2] 



(n - k)(k + n+ 1) 



Pc(*) 



E 

fc=0 



2n - 4fc - 1 
(n-fc)(2fc+l) 



P>- 



2fc-l 



(z). 



(3.46) 



(3.47) 



(3.48) 



A glance at Eq. (3.48) reveals that W n -\{z) is the well-known Christoffel's polynomial (cf, e.g., [14, 
page 153]), in terms of which one has [23,25] 



Qn{z) = \Pn(z) lni±| - W„_i(*). 



(3.49) 



On combining Eqs. (3.45), (3.46), and (3.49), we obtain the sought relationship 3 [cf Eq. (1.1)] 

ap<t(z) 



= iP(n + l)P n (z) + Q n (z). 



(3.50) 



The same result may be obtained if one couples Eq. (3.2) with Eq. (4.1), both particularized to 
the case of m — 0. 



3 The counterpart relationship (see [23, page 178]) 

dQZ(z) 



~iw + iP(n + l)]Q n (z) 



M=0 



follows straightforwardly from Eq. (2.30) after differentiating the latter with respect to /j, and setting then fj, = 0. 
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3.2.3 The case of to > n 

In this case, the function P™(z) vanishes identically [cf Eq. (2.10)], so that Eq. (3.4) becomes 

dP%(z 



= U™(z) (m>n). 



fi—m 



With the help of Eq. (3.37), from Eq. (3.6) we find 



(-)'■ 



z-l 



m/2 n 

£(- 

fe=0 



k (k + n)\(m - k - 1)! f z-l 
' k\(n-k)\ 



Comparison of Eq. (3.52) with Eq. (2.20) reveals the relationship 4 

U™(z) = (-) rn (n + to)!(to - n - l)\p- m (-z) (m > n) 
Consequently, from Eqs. (3.51) and (3.53) one obtains 

dP%(z) 



(3.51) 



(m > n). 



(-) m (n + to)!(to - n - l)\p- rn (-z) (to > n). 



(3.52) 



(3.53) 



(3.54) 



3.3 Evaluation of [dP%{z)/d^-_ 
3.3.1 The case of < to < n 

For ii = —m, Eq. (3.4) becomes 



dP£(z) 



H=—m 



= lp- m (z)ln^±± + U- m (z) 



Since it holds that 



dPX(z) 



dP-»(z) 



d/j, 



(0 ^ to ^ n). 



from Eqs. (3.55), (3.2), and (3.35) it may be inferred that 

U~ m (z) = (-) n+1 [ n ~ ^ U?(-z) + Mn + to + 1) + 4,{n - to + l)]P» m (. 
(n + my. 

(0 < to < n). 



(3.55) 



(3.56) 



(3.57) 



One may couple this formula with Eqs. (3.38) to (3.43) and exploit, whenever necessary, Eqs. (2.11) 
and (2.12) to obtain various explicit representations of U~ m {z) with < m ^ n. 



3.3.2 The case of to > n 

From Eq. (3.4) we find 

dP£{z) 



fj,=—m 



= Ip-™ ( *) In J±l + [/-">(*) (to> 



n . 



(3.58) 



4 The result (3.53) may be also deduced from Eq. (3.34). Indeed, it is seen that in the case considered here the 
summand in the latter equation vanishes identically and only the first term on the right-hand side survives. Further, 
in virtue of Eqs. (2.9) and (3.37) (with k replaced by n), we have 

lim Mn -n+ l)P£(z) = (-) m (n + m)!(m -n- l)\p- m (-z) (m > n), 

— >m 

hence, Eq. (3.53) follows. Still another way of arriving at Eq. (3.53) is to use Eq. (3.5). 
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Equations (3.6), (3.10), (3.12), and (3.13) supply us with the following alternative representations 
of [/-•"(*): 



z + 



fe=0 



fc!(fc + m)!(n-fc)! 



t/- m (z) - [^(n + m+l) + ^(m-n)]p- m (z) 



1 (z - 1 

(n + m)!(m-n - 1)! + 1 



m/2 



E 

fc=0 



(fc + n)\(m - k - l)\ip(m - fc) ( z + 1 
fc!(n-fc)! 



(to > n), 



(3.60) 



C/- m (z) - ^(m-n)p- m (z) 



z - 1 



><E(- 

fe=0 

(to > n), 



(m — n — 1)! \ z + 1 
fc (fc + m — n — 1)! 



i/2 



Z+l 



fc!(fc + m)!(n-fc)! 



[tp(k + m + 1) — ip(k + to — n) 



z-1 
z + l 



(3.61) 



U~ m (z) = 1 p(m-n)p- m (z) + (-r 

n 



z-1 

(to — n — 1)! ^z + 1 
(m-fc-1)! 



t/2 



Z - 1 



[V>(n + m — fc + 1) — -</)(m — fc)] 



fc=0 

(m > n). 



k\{n — fc)!(n + to — fc)! 
Moreover, with the aid of Eq. (2.23), from Eq. (3.33) we obtain 

(fc + to)!(to - fc - 1)! 



z + l 
z-1 



(3.62) 



U- m {z) = 4>{n + m+l)p- m {z) 

n-1 

+E(-) fe+ "— 

(to > n). 



2fc + 1 
(n - fc)(fc + n+ 1) 



(n + to)!(to — n — 1)\ 



(3.63) 



3.4 The derivative dP£(x)/dpi 



Since [cf Eq. (2.24)] the definition of P£(x) with -1 < a; < 1 differs from that of P£(z) with 
z G C \ [—1, 1], the derivative dP^(x)/d/i requires to be considered separately. 

Differentiating Eq. (2.24) with respect to fi, making use of Eq. (3.35) and exploiting the fact 
that 

.x + l±i0 = l + .x a:-l±i0 = c ±i7r (l-x) (-1 < or < 1) (3.64) 



results in 



where, in analogy with Eq. (2.24), we define 



K(x) = - 



^ /2 U%{x + iO) + c-'™»l 2 U£{x - iO) 



3 ±i7T(U 



/ 2 [/£(x±i0). 



(3.65) 



(3.66) 



Various representations of Ufi(x) and U^ m (x) arise if one couples Eq. (3.66) with the formulas for 
U%(z) and U^ m (z) derived in sections 3.1 to 3.3 (now particularized to the case of z = x ± iO). 
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4 Applications 



4.1 Construction of the associated Legendre function of the second kind 
of integer degree and order 

As the first application of the results obtained in the preceding section, below we shall construct 
several explicit representations of the associated Legendre function of the second kind of integer 
degree and order. 

4.1.1 The functions Q± m {z) in the case of < m < n 

If in Eq. (2.26) we let \i approach ±m, after applying the l'Hospital rule to the right-hand side, we 
obtain 



Qt m {*) 



1 dPX(z) 



2 dfi 



(-)" dP£{-z) 



(0 ^ m ^ n). 



(4.1) 



fi—^zm 



fi—^zm 

Combining this with Eq. (3.3) gives the following Rodrigues-type representation of Qn m ( z ) with 
< m ^ n: 



+ 



+ 



z-l 



2 n + 1 (n + m)\ + 1 



z + 1 



m/2 



1 



2«+i( n =Fm)! \z - 1 
(0 ^ m ^ n). 



2/2 



dz™ 

d" 
dz™ 



(z- l) n - m (z + l) n+m ln 
(z-l)" +m (z + l)"- m ln 



z + 1 

z-l 

z + 1 
z-l 



For m = 0, Eq. (4.2) reduces to the well-known (e.g., [25, Eq. (8.836.1)]) formula 

z + 1 



Q ^=-l p ^z^i + ^£ 



(z z -l) n \n 



z - 1 



Upon making use of Eqs. (3.4) and (2.11), Eq. (4.1) implies that 



1 -y±m, x, z + 1 



with 



Qt m (z) = ~ 2 P'i W^[(z) (0 ^ m ^ n), 
W^[(z) = \ [(-) n Ut m (-z) Ut m (z)] (0 < m < n). 



From Eq. (4.5) it is seen that 

W£!i(-z) = (-T +1 Wt+\{z) (0 < m < n), 



(4.2) 
(4.3) 

(4.4) 
(4.5) 

(4.6) 



i.e., the parities of W n ^[(z) are opposite to those of P^' m (z). In view of the symmetry relation 



(4.7) 



which is the consequence of Eqs. (4.5), (3.57), and (2.11), henceforth we shall discuss W^Li(z) 
only. 

A number of explicit representations of W™_ 1 (z) may be obtained if one couples Eq. (4.5) (with 
the plus signs chosen) with Eqs. (3.38) to (3.43). For instance, if U™(—z) is calculated using Eq. 
(3.39), while U™(z) is taken in the form (3.38), then Eq. (2.32) is recovered. Conversely, if U™(—z) 
is obtained from Eq. (3.38) and U™(z) from Eq. (3.39), this results in Eq. (2.33). If in Eq. (4.5) 
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both U™(—z) and U" L (z) are obtained from Eq. (3.38), one recovers the Brown's formula (2.34), 
while if both U™(—z) and U™(z) are derived from Eq. (3.39), one arrives at 



+ 



1 (n + m)! (z- l\ m/2 ^ (fc + n)!^(fc + m + l) (z-l 



2{n-m)\ \z + l 



(-)" (n + m)! /z + 1 



^ k\{k + m)\{n - k)\ 



k=0 
m/2 



2 (n — m) 
(0 < m ^ n). 



^ fe (k + n)\ip(k + m + 1) + 1 
/c!(fc + m)!(n-/c)! 



(4.8) 



That of the two Snow's representations which corresponds to the choice of the upper signs in Eq. 
(2.35) follows if in Eq. (4.5) U™(-z) is found from Eq. (3.41) and U™(z) from Eq. (3.40); the 
other one is obtained if the roles played by Eqs. (3.40) and (3.41) are interchanged. 

An apparently new representation of W™_ 1 (z) follows if in Eq. (4.5) both U™(—z) and U™(z) 
are found from Eq. (3.33). Proceeding in this way, with the aid of Eq. (2.11), one finds 



1 (n + m)! 

2 (n-m)! 

— rn- 

E 



m—l 



2fc + 1 



k=0 

!-(-)* 



(n- fc)(fc + n+ 1) 

2fc + 2m + l 



k=0 



(n — m — k)(k + n + m + 1) 



i-H 



k+n+m 



k\(n + m)! 



(fc + 2m)!(n — m)\ 



Pr +m {z) (0<m^ n). (4.9) 



It is seen that the summand in the second sum on the right-hand side of the above equation 
vanishes if k + n + m is even and therefore Eq. (4.9) may be rewritten more suitably as 



wm / \ 1 {n + m)\ 

" 2(^1)! 



m—l 



2A:+ 1 



fc=0 



(n- fc)(fc + n+ 1) 



-[P fc - m (-z)-(-)"Pr"(z)] 



+ E 



l-(-) 



fc+n+m 



2fc + 2m + l 



fc=0 



1 + 



2 (n — m — k)(k + n + m + 1) 

k\(n + m)! 



(k + 2m)!(n - m)! 
Manipulating with the second sum, Eq. (4.10) may be cast into 

m—l 



Pk+mW (0^m<n). 



(4.10) 



W^iz) 



1 (n + m)! fc 2fc+l 

2(n-m)!^ J (n - fc)(fc + n + 1) 



+ 



int[(n-m-l)/2] 

1 >— ^ 2n — 4fc — 1 

? ^- (n-fc)(2fc + l) 



fc=0 



1 + 



(n + m)!(n — m — 2k — 1)! 



(n-m)!(n + TO-2fc-l)! 



(0^m<n). (4.11) 



It is seen that in the case of m = Eqs. (4.10) and (4.11) reduce to Eqs. (3.47) and (3.48), 
respectively. 

4.1.2 The functions Q™(z) and Q™ n _i(z) in the case of m > n 

Proceeding as in section 4.1.1, from Eq. (2.26) we obtain 

(-)" dP£{-z) 



1 BP^(z) 

Az) = 2^r 



fi—ra 



(m > n). 



(4.12) 



fj.—m 
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If use is made here of Eq. (3.54), this yields Eq. (2.37). Next, for the function Q™ n _ 1 (z), from Eq. 
(2.27), with the aid of the l'Hospital rule, we have 



Q™ n - 1 (z) = -Q: i ( Z ) + 



dP%{z) 



dp 



(ra > n). 



(4.13) 



From this, after merging with Eqs. (2.37) and (3.54), we find 

Q m n -i(z) = ^-(n + m)!(m - n - 1)! [p- m (-z) + (-) n P- m {z)] (m > n). (4.14) 



4.1.3 The functions Q± m (x) and Q™„_i(a;) 

Inserting Eq. (4.4) into Eq. (2.38), the latter being particularized to fi — ±m, with the help of Eq. 
(2.24) we obtain 



QTW = In — W^{x) (0 < m < n), 



(4.15) 



where we define 



W^(*) = 



(-)' 



gTim/J^^ + iQ ) + e ±inrn/2 w ±rn( x _ i0 )l ( <; TO < n) 



(4.16) 



It follows from Eqs. (4.16) and (4.6) that the functions W^H{{x) possess the reflection property 

W^(-x) = (-) n+m+1 W±™(x) (0<m<n), (4.17) 
while if Eq. (4.16) is coupled with Eq. (4.7), this results in 



W~^(x) = (-r{^^WC-iM (0 < m < n). 



(4.18) 



The latter relationship allows one to focus on W^ l _ 1 (x) only. 

Exploiting the expressions for W™_i(z) found in section 4.1.1, with no difficulty one may con- 
struct counterpart formulas for W^ l _ 1 (x). For instance, if Eq. (4.16) is combined with Eq. (4.11) 
and use is made of Eq. (2.11), this yields 



m — 1 



2k + 1 



2(n-m)!^ v (n-k)(k + n+l) 



+ - 



1 int [( ^-l)/2] 2n _ Ak _ l 



E 



k=0 



(n- fc)(2fc + 1) 



1 + 



(n + m)!(n — m — 2fc — 1)! 



n-2fc- 



(n-m)!(n + m-2fc- 1)! 
If Eq. (2.37) is used in Eq. (2.38), one finds 

On («) - + m)\(m - n - 1)! [P- ro (-z) - (-)" +m p- m (x) 

while if Eq. (4.14) is coupled with Eq. (2.39), this gives 



^(x) (OsCmsCn). (4.19) 



(m > n) 



(4.20) 



(-)' 



■(n + my.(m-n-l)\[p- m (-x) + (-) n+m p- m (x)] (m > n). (4.21) 
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4.2 Evaluation of [d 2 P^{z)/diJi 2 ] lJL=m for m> n 

In this section, we shall show that if to > n, then the knowledge of [dP%{— z)/dfi] ll= - m allows one 
to find [d 2 P%{z)/diiX =m . 



We begin with the observation that from the identity 

V>(0 = ^(i - - cos(7rC)r(c)r(i - c), 



(4.22) 



which may be easily deduced from elementary properties of the digamma and gamma functions 
[14,20,23,25,30], one finds that 



ip(n - n + 1) = V(m -n) + (-)" cos(n fj,)T(fj, - n)T(n - fj, + 1). 
The identity (4.23) allows one to rewrite the relation (3.2) in the form 
dP-^z) 



(4.23) 



.^(n + fx+l) + ^-n)}p-^z) 



+ (-) 



n T(n-fi + l) 
r(n + /i+l) 



cos(w^)T(fi - n)T(n + p. + l)P~»{z) 



(4.24) 



hence, it follows that 

dP-^(z) 



fi=m 



■[ip(n + to + 1) + ip(m - n)]P~ m {z) 
(-)' 



(n + to)! iu— >m 
~dP%{-z) 



lim iT(n-fj, + l) 

- cos{n(i)T(fj, - ri)T(n + fj, + l)P~^{z) 



(to > n). 



(4.25) 



The limit in the above equation may be evaluated with the aid of the l'Hospital rule [for this 
purpose, the property (3.37) proves helpful] and one finds 



dP-»{z) 



djj, 



+ 



2[iP(n + m + 1) + V(m - n)]p- m {z) 



dP-^(z) 



(n + to)!(to — n — 1)! d[i 2 



(m > n). (4.26) 



If Eq. (4.26) is solved for [d 2 P£(— z)/dfi 2 ]^ =m , then use is made of Eq. (3.56) and z is replaced by 
—z, one eventually arrives at 



d 2 P£(z) 



fi—m 



(-) m 2(n + TO)!(TO-n-l)! 



x \ [tp{n + to + 1) + yj(m - n)]P" m (-z) - 



dP£(-z) 



fi— — ra ^ 



(m > n). 



(4.27) 



Several explicit representations of [d 2 P^(z)/dfj 2 ] fl=m with to > n, not listed here, follow if one 
combines Eq. (4.27) with the results of section 3.3.2. 
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n 



4.3 Evaluation of [dQ^(z)/dfj] lx=m and [dQt n -i(z)/dfj] fi=m for m > 

The results of section 4.2 may be exploited to express the derivatives [dQ^(z)/dfi\ ll — m and 
[dQ-n-iW/dA^m for m > n in terms of [dP£(±z)/dfi\ lt =- m . 

To show this for [dQ^(z)/dfi] li - m , we differentiate Eq. (2.26) with respect to (i, obtaining 



djj, 



sin(nfi) 



1 ^ dP£{z) _ (-)" j^ dPXj-z) 



.coB{w)Q»{z) + -j«»-j±J--^ 



dfx 



(4.28) 

In the limit \i — > to > n, with the help of the l'Hospital rule and of Eq. (4.12), from Eq. (4.28) we 
obtain 



dQZ(z) 



(m > n), 



5Q£(z) 



1 &P>;(z) 



fi—m 



2 



(j,=m 



(-)" 2 P£(-z) 
2 <V 



fj,=m 



(4.29) 



hence, it follows that 



9yU 



1 d 2 P' J -(z) 

= m Qn \ Z ) + 7 



/j,=m 



in v ' 4 a M 2 



(-) n ^(- 2 ) 



ctyt 2 



(to > n). 



(4.30) 

If the last two terms on the right-hand side of Eq. (4.30) are transformed with the aid of Eq. (4.27) 
and use is made of Eq. (2.37), this yields the sought relationship 



dQZ(z) 



fi—rn 



[in + ip(n + to + 1) + ip(m - n)]Q™{z) 



{-)' 



■(n + m)\{m — n — 1)! 



dP%{-z) 



8fM 



fj,=—m 



l_i— — m 



(4.31) 



2 

(to > n). 

If one starts with Eq. (2.29) rather than with Eq. (2.26), after movements almost identical to 
these presented above one arrives at 



= [in + -4>{n + to + 1) + ip(m - n^Q^^z) 



(-)' 



2 

(to > n). 



(n + m)!(m — n — 1)! 



dP£(-z) 



+ H 



ndP^(z) 



l_i— — m 



f_t— — m 

(4.32) 
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